Active polymer translocation through flickering pores 
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Single file translocation ol a homopolymer through an active channel under the presence of a 
driving force is studied using Langevin dynamics simulation. It is shown that a channel with sticky 
walls and oscillating width could lead to significantly more efficient translocation as compared to 
a static channel that has a width equal to the mean width of the oscillating pore. The gain in 
translocation exhibits a strong dependence on the stickiness of the pore, which could allow the 
polymer translocation process to be highly selective. 
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Introduction. — The translocation of a polymer through 
a pore is important in the context of many biological pro- 
cesses such as the transport of RNA through a nuclear 
membrane pore [T] and the injection of viruses. Its var- 
ious technological applications such as drug delivery [2], 
rapid DNA sequencing [2H1] and gene therapy has led to 
several recent experimental [JHZj and theoretical studies 
[SIBU]. Most theoretical work in this field has focussed 
on the underlying physics of the translocation process 
O El [H] and the effects of the pore-polymer interactions 
[Tot [Mj [T8] , structure of the pore [17] , crowding [19] and 
confinement effects on the dynamics of translocation. Ex- 
perimental studies on electric field driven translocation 
of DNA and RNA molecules across a-hemolysin channels 
[1] prompted the introduction of additional driving forces 
to aid translocation. In most of these studies the pore 
is considered static with the polymer always experienc- 
ing a constant confinement during its translocation from 
the CIS to the trans side. However, there are a number 
of biological examples, such as the twin-pore translocase 
complex in the inner membrane of mitochondria |21) and 
the nuclear pore complex |22j , where it is known that the 
width of the channel effectively changes during the course 
of translocation. Inspired by these examples, we set out 
to study the generic effect of such temporal modulations 
on the efficiency of polymer translocation, using a simple 
coarse-grained model. 

We find that the translocation of a polymer through a 
narrow channel with a width that oscillates with a given 
frequency can be significantly enhanced when compared 
to a static pore. The time of translocation is sensitive to 
the initial condition, and the driving force and stickiness 
of the channel walls significantly affect the translocation 
at the limit of high frequencies. 

Model. — In our simulation, we model the polymer as 
a bead-spring chain. The polymer beads experience 
an excluded volume interaction modeled by a repulsive 



Lennard- Jones (LJ) potential of the form U^^n^ir) 
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strength of the LJ potential. The monomers of the 
chain experience an additional interaction modeled as 
a finitely extensible nonlinear elastic (FENE) spring: 
U^™^{r) = -ifci?2ln(^l _ where k is the spring 
constant and R is the maximum separation between con- 
secutive monomers along the chain. We consider a two- 
dimensional (2D) geometry where the pore is modeled as 
being made up of monomers of diameter a as shown in 
Fig. [l^ (inset). The walls perpendicular to the pore 
are taken to be long enough to avoid the crossing of 
the polymer ends. The interaction of the polymer with 
the wall, U^ii, is modeled by the same excluded vol- 
ume interaction that exists between the polymer beads; 
thus, = C/j^im. The interaction of the polymer 

with the pore is modeled by the Lennard- Jones poten- 



tial U}^i{r) = 4ep 



where cr is the diameter of a bead and e gives the 



'pmvw - -^pm Vr/ \rJ ^ - ^-Sct and 

for r > 2.5(7. Inside the pore, the polymer experiences 
an external driving force F^xt — F^i directed along the 
pore axis. A Langevin dynamics algorithm is used to 
integrate the equation of motion of the polymer beads 
mVi — —^Ui + Foxt — C'^i + flii where m is the monomer 
mass, C/, = C/^i, + U^h^^ + C^wm + U^i is the total 
potential experienced by a bead, C, is the friction coef- 
ficient, Vi is the monomer velocity, and rjj is the ran- 
dom force satisfying the fluctuation-dissipation theorem 
{rl^{t) ■ Vj{to)) = AkBTCS,jS{t - to). 

In our model, e, cr, and m set the units of energy, 
length, and mass, respectively, which result in a unit of 
time as (mtr^/e)^/^. Using these units, the dimensionless 
parameters oi R — 2, k — 7, k-QT — 1.2, and C, = 0.7 have 
been chosen for the simulations, in accordance with ear- 
lier simulation studies of polymer translocation dynamics 
[m [15] . The length of the pore is fixed at L = 5, and the 
length of the polymer is TV = 32 unless otherwise speci- 
fied. We consider the case where the width of the pore, 
W{t), is allowed to oscillate harmonically, with frequency 
uj — 21: 1 T and amplitude Wa^ about an average width 
Wq = (I^(i))T' namely W(t) = Wq VTa sin (wt (/i) , 
where is an initial phase. Note that Wmin = W^o — 
and Wijiax = W^o + ^ A are the minimal and maximal 
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FIG. 1: (a) Gain as a function of the scaled frequency, for F — 1.0. Top inset is a schematic of the simulation geometry. 
Bottom inset shows an oscillation cycle of the pore, (b) High frequency limit of the gain for different values of tpm and F. The 
points are simulation data for oscillating widths and the lines are extracted from an instantaneous static width approximation 
(see text). Top inset shows the translocation velocity as a function of the static width. Bottom inset shows the velocity as a 
function of time during the oscillation cycle, (c) The translocation time as a function of width for static pores, for = 1.5 
and F = 1.5. The points are simulation data. Inset shows the same plot for different forces. 



widths of the pore, respectively. Due to the strongly 
repulsive excluded volume interaction between beads, re- 
ducing the pore width below a certain value could re- 
sult in the breaking of the bonds between neighboring 
monomers in the polymer. We choose Wmin such that in 
a static pore with this width a polymer would be trapped 
without breaking up. Moreover, the requirement of single 
file translocation of the polymer (i.e. no hairpin struc- 
tures), would limit W^max- In accordance with these re- 
strictions, wc choose Wq — 2.25 and Wa — 0.25, so that 
the pore oscillates between Wmin = 2.0 and Wmax = 2.5. 
(The value for the relative change in width might ap- 
pear too large if regarded as a conformational change in 
a protein. However, this is inflicted by the particular 
model of rigid spheres used in our coarse-graining, and 
could be much smaller for more realistic descriptions.) 
We have checked that our results are not sensitive to the 
particular value of the amplitude, so long as it satisfies 
the requirements that it closes the channel at minimum 
width and allows relatively free passage of the polymer 
at maximum width. The driving force inside the pore 
varies between 0.5 and 2. Initially, the first monomer 
of the polymer is held fixed at the entrance of the pore 
while the other beads are allowed to fluctuate. After al- 
lowing sufficient time for the polymer configuration to 
equilibrate, the first monomer is released and the time 
that elapses between the entrance of this monomer into 
the pore and the exit of the last monomer is measured. 
This gives the translocation time of the polymer through 
the pore. The time step in our simulations is chosen as 
At = 0.01 and the averaging is done over 2000 successful 
translocation events. 

Results. — We examine the efficiency of the transloca- 
tion process by comparing the average time of transloca- 



tion for the oscillating pore Tqsc with the average translo- 
cation time for a static pore (of width Wq) Tgtat- We 
first focus on the (f> — case. In Fig. [T^, the gain in 
translocation rate defined as 77 = Tstat/Tosc is plotted as 
a function of the dimensionless frequency w = Tstat/T. 
We find that the translocation rate is enhanced for the 
oscillating pore as compared to the static pore. For the 
iV = 32 case, there are two distinct peaks at a) ~ 0.5 
and 1.5, with the translocation time at w ~ 0.5 almost 
half of that for the static pore. This behavior can be 
understood by noting that during the first half period 
of oscillation of the pore (between t = and t = T/2), 
the width of the pore is always greater than the aver- 
age width, 14^0, which is the width of the static pore. 
When the oscillation frequency is sufficiently small, then 
the polymer only experiences this half of the cycle before 
the completion of the translocation process. The average 
translocation time therefore increases steadily for small 
frequencies. Beyond the critical oscillation frequency cor- 
responding to uj 0.5, the polymer experiences the effect 
of the second half of the cycle where the pore width is 
always smaller than the static pore width, Wq. The aver- 
age translocation time increases, thus lowering the gain 
substantially, until it reaches a minimum corresponding 
to w ^ 1. For higher oscillation frequencies, the polymer 
starts experiencing the next half of the cycle between 
t — T and t = 3T/2, when the pore width is again larger 
than the static pore. Therefore, we observe a second peak 
in the gain corresponding to Tgtat ~ 5T/2. This peak is 
much smaller in magnitude because at higher frequencies 
a substantial amount of time is spent by the polymer in a 
trapped state. As the frequency is increased further, the 
gain becomes insensitive to the frequency and develops 
a plateau [23]. Note that as the frequency tends to zero. 
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FIG. 2; Translocation time distribution (top) as compared 
with the distribution for the static pore with average width 
(soUd hne), and average number of translocated monomers as 
a function of time (bottom) at a; = 5.19, for F = 1.0 and 
enm = 1.0. 



the pore is essentially static and the polymer transloca- 
tion time reduces to the translocation time for the static 
pore and the gain approaches unity. We have observed 
that the behavior of the system is robust as the length of 
the polymer is changed. For example, Fig. [T^ shows that 
for N = 128 all the features are identical to the = 32 
case, except that two additional higher frequency peaks 
are observed. 

We find that the gain increases as epm (that repre- 
sents the stickiness of the pore) is increased, although 
the general features — the primary and secondary peaks 
at intermediate frequencies and the saturation at high 
frequencies — are preserved (see Fig. [T^). While the low 
frequency limit of the gain is universal, we observe that 
the saturation value of the gain at high frequencies de- 
pends on the external force and the stickiness of the pore, 
as shown in Fig. ^p. 

Since the translocation is controlled by the constraint 
set by the time varying width of the pore, studying how 
the translocation time for a static pore depends on the 
width could help us understand the observed behavior 
[Mj. In Fig. [ij;, the average time for the translocation of 
a polymer through a static pore is plotted as a function of 
the pore width. One can identify two distinct behaviors, 
namely, a very sharp decrease for Wmin = 2.Q < W < 
2.25, which crosses over to a regime with relatively slower 
decay for 2.25 < W < L = 5.0. We find that the former 
regime is controlled by the time the polymer takes to 
transfer through the pore, whereas the latter is controlled 
by the time it takes to escape the pore. The inset in Fig. 
[l]: shows the same behavior for different values of the 
external force. 

One can gain further insight by using a scaling argu- 



ment for a confined polymer within the blob picture [53] . 
A confined polymer that fills the entire pore of length 
L breaks up into L/W blobs of uniform size W. How- 
ever, during escape the polymer is split into a fraction 
that fills a part of the pore of length L — y, leaving the 
part of length y empty, and a part that has gone out- 
side. The entropic penalty for the partial confinement 
of the polymer inside the channel is ~ A:bT(L — y)/W . 
The monomers that are inside the channel will experience 
an external force F. Therefore, displacing the polymer 
by a length y inside the channel will result in a gain 



of -F{L^ 



y 



2Ly)/2W in mechanical energy. The 



interaction between the monomers and the sticky walls 
of the pore should also be taken into account. Since 
the LJ attraction is short ranged, we can estimate this 
contribution by assigning an energy gain of —epm to ev- 
ery monomer that is in the vicinity of the channel walls. 
Counting these monomers, we find this "adsorption en- 
ergy" as -2ep^[{L - y)/W]W'^/^{l/W). Therefore, the 
total free energy for escape (up to a constant) is given 



as T{y) = 



FL 
W 



2W 



where ci 



and C2 are constants of order unity. We can now treat 
the escape problem as diffusion across a one-dimensional 
effective potential barrier, and calculate the mean first 
passage time [TU]. The resulting escape time would lead 
to a similar trend as in Fig. [l]:. 

We can now calculate the translocation velocity, de- 
fined a,s V = N/Tstat, for a static pore as a function of 
the width, and use the time dependence of the width to 
extract the translocation velocity as a function of time 
over a full period. These are shown as insets in Fig. [i]d. 
Integrating the velocity over a time range that is just 
long enough to allow for the translocation of the entire 
polymer, we get the translocation time for an oscillating 
pore within this approximate scheme. The solid lines in 
Fig. [T]3 show the result of this calculation for the high 
frequency saturation value of the gain for different val- 
ues of F and ep,„. The instantaneous static pore picture 
thus seems to provide a reasonable account of the phe- 
nomenon, although it systematically underestimates the 
gain due to the absence of noise. The asymmetry ob- 
served in the two half cycles of the translocation velocity 
is a direct consequence of choosing the average width Wo 
to be the crossover point between two different regimes 
in the translocation time versus width plot in Fig. 

It is instructive to examine the distribution of the 
translocation times in the high frequency regime. As 
shown in Fig. [2j we find that the distribution consists 
of a series of peaks separated by distinct minima. These 
minima correspond to the points in the oscillation cycle 
when the width of the pore is at its minimum, Wmin- 
We can further explore the detailed dynamics of polymer 
translocation by looking at the number of translocated 
monomers, iVtrans, i-e- the number of monomers that 
have left the pore as a function of time. The plot shown 
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FIG. 3: Gain as a function of the scaled frequency for dif- 
ferent values of cf>. The plots correspond to F = 1.0 and 
Epm = 1.0. Inset shows the corresponding width oscillation 
cycles. 

in Fig. [2] is calculated by averaging over all successful 
translocation events. iVtrans as a function of time alter- 
nates between flat and linearly increasing regions. The 
flat regions correspond to the time periods in the oscil- 
lation cycle when the pore width varies between Wmin 
and Wo, when the polymer has very little space to ma- 
noeuvre itself and is essentially trapped. Therefore, the 
monomers are not able to leave the pore, and hence, 
-^trans docs not change. During the periods in the oscilla- 
tion cycle when the width varies between and Wmax, 
the polymer is no longer trapped and the monomers can 
escape the pore resulting in an increase in iVtrans- Note 
that an increase in the frequency of oscillation uj causes 
an increase in the number of peaks in the translocation 
time distribution. 

In our analysis so far we have assumed that the pore 
width at the beginning of the translocation process is 
the mean width of the oscillation cycle. Figure |3] shows 
the effect of the initial phase (j) on the gain as a func- 
tion of the oscillation frequency. We observe that at low 
frequencies the initial phase strongly affects the translo- 
cation gain by controlling the effective width during the 
translocation period, in agreement with the picture de- 
scribed above. However, as the initial phase only affects 
the first couple of cycles, we expect that at sufficiently 
high frequencies the translocation process becomes in- 
sensitive to it. This is confirmed by the observation that 
the saturation value of the gain at high frequencies is in- 
dependent of the initial phase, as Fig. [3] demonstrates. 
The robust high frequency behavior of the system has 
remarkable implications. It suggest that even stochas- 
tic fiickering of the pore will result in a steady average 
translocation with a gain that can be tuned via F and 



Epm, as long as T <C Tstat, where T is the characteristic 
period of the random opening and closing of the pore. 
Moreover, we can infer from Fig. [ifj that the translo- 
cation rate for such random flickering pores is extremely 
sensitive to the stickiness of the pore, which could bring 
about the possibility of a high degree of robust selectivity 
in the translocation process. 

It will be interesting to probe whether such poten- 
tial for selectivity is already exploited in biological active 
pores. Other interesting factors could be the flexibility of 
the polymer and the strength of the noise, which could 
potentially lead to stochastic resonance of the polymer 
|26) . The design rules that can be obtained from our 
study could also be used in fabricating highly selective 
synthetic active pores, which might be able to perform 
tasks such as sequencing. 
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